A model of hyphal tip growth involving microtubule-based transport 
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We propose a simple model for mass transport within a fungal hypha and its subsequent growth. 
Inspired by the role of microtubule-transported vesicles, we embody the internal dynamics of mass 
inside a hypha with mutually excluding particles progressing stochastically along a growing one- 
dimensional lattice. The connection between long range transport of materials for growth, and the 
resulting extension of the hyphal tip has not previously been addressed in the modelling literature. 
We derive and analyse mean-field equations for the model and present a phase diagram of its steady 
state behaviour, which we compare to simulations. We discuss our results in the context of the 
filamentous fungus, Neurospora crassa. 

PACS numbers: 87.10.-(-e, 87.16.Ac, 87.16.Ka , 05.40.-a 



Biologically, fungi are distinct from both plants and 
animals. In addition to their intrinsic interest, they im- 
pact immensely on human affairs and on the ecosystem 

Key to the evolutionary success of fungi, is their unique 
mode of growth. Filamentous fungi grow by the polarised 
extension of thread-like hyphae, which make up the body, 
or mycelium, of a fungus. Except for branching (which 
initiates new hyphae) the site of growth is localised to a 
single region at the tip of each elongating hypha. 

There are many theoretical models for the growth of 
fungal colonies and of single hyphae (reviewed in [1, 
Most models of single hypha growth concentrate on bio- 
mechanics 0, S|- Of more interest for us here is the 
"vesicle supply centre" (VSC) model [1,0] j in which raw 
materials for growth are packaged in secretory vesicles 
and distributed to the hyphal surface from a single "sup- 
ply centre" (often identified with an organelle complex 
known as the Spitzenkorper, or apical body 0, 0|) situ- 
ated within the growing tip. This model is capable of pre- 
dicting the shape of hyphal tips; but the speed of growth 
(equivalent to the speed of the VSC) is an input parame- 
ter. Moreover, all transport processes are subsumed into 
a single rate of vesicle supply at the VSC. We are aware 
of just one model that takes explicit account of trans- 
port along the growing hypha [10]. A major interest of 
this early work, however, was the initiation of branching; 
these authors did not relate vesicle transport to growth 
velocity. This latter issue remains poorly understood. 

In this work, we propose a simple one-dimensional 
model which makes an explicit connection between the 
long-distance transport of building materials along a hy- 
pha and the resulting extension as they are delivered to 
its apical site of growth. 

It is a highly idealised model, encompassing many 
complicated biological processes (many of which are still 
poorly understood) with two key parameters: the rate at 
which vesicles enter the system and the efficiency with 



which they extend the length of the hypha. We demon- 
strate that by altering these rates, steady states can be 
attained whereby the hypha is extending at a constant 
speed while being supplied with materials far behind the 
tip. Our model has features in common with 10]. Like 
[ig ]. we use computer simulations, but unlike we can 
also make analytical progress because of recent advances 
in non-equilibrium statistical physics. 

Our model is inspired by the well-known Totally Asym- 
metric Simple Exclusion Process (TASEP). This is a 
one-dimensional lattice along which particles progress 
through stochastic directed motion. No more than one 
particle may occupy each lattice site at any given time. 
The TASEP was originally introduced as a lattice model 
of ribosome motion along mRNA llj and recent variants 
have been widely used to model the collective dynamics 
of molecular motors 12, 13, 14, 1^. The application of 



this and other classes of statistical mechanical models to 
many kinds of "biological traffic" has recently been re- 
viewed [l^ . The TASEP is also widely studied in its own 
right as a fundamental model of non-equilibrium statisti- 
cal mechanics 17[ , in particular as a simple driven diffu- 
sive system exhibiting non-equilibrium phase transitions 
[3] between different macroscopic density and current 
regimes [lo!, ^2d\. Our work contributes to the study of 
both "biological traffic" and non-equilibrium phase tran- 
sitions. 

In constructing the model, we introduce a new feature 
into the TASEP: particles reaching the end of the lat- 
tice act to extend it. We ask whether a constant input 
rate far from the growing end of the lattice can gener- 
ate steady state lattice growth and if so, how the growth 
velocity depends upon the system parameters. We find 
that, as in the TASEP, different macroscopic regimes ex- 
ist in the growing system, with different forms for the 
growth velocity, and non-equilibrium phase transitions 
between these regimes. We ask whether these steady 
states may be related to the growth states observed in 




FIG. 1: A simplified diagram of the process of polarized se- 
cretion resulting in tip growth in a fungal hypha. Secretory 
vesicles are transported by kinesin molecular motors (not to 
scale) along microtubules towards their growing "plus" ends. 
The vesicles accumulate within the so-called Spitzenkoi^per 
at the hyphal tip before fusing with the apical plasma mem- 
brane. The secretory vesicles deliver membrane proteins and 
lipids, cell wall synthesizing enzymes and possibly cell wall 
precursors to the growing fungal tip. Hyphal widths typically 
vary between 5 and 15 nm. 

fungal hyphae. 

To arrive at the model we appeal to the popular be- 
lief that the long distance transport mechanism within 
a hypha is provided by kinesin molecular motors, "walk- 
ing" along microtubule filaments which run length-ways 
through the hypha Once at the tip, vesicles fuse 

with the plasma membrane, resulting in localised exten- 
sion of the hyphal wall [1, 0, [llj . Fig. [T] summarises this 
process diagrammatically. 

A kinesin motor with cargo progressing toward the tip 
will likely attach and detach between a number of micro- 
tubules on its way, and is then thought to be transferred 
onto actin filaments at the VSC, before being distributed 
to the hyphal tip In our simple model however, we 
shall bypass the details of the vesicle trajectory and con- 
centrate instead on a coarse-grained description repre- 
senting what is essentially a continuous directed move- 
ment of vesicles toward the hyphal tip. We thus assume 
here simply that a continuous driving force is applied to 
the vesicles, all the way to the point of fusion. 

We construct the model by supposing that a hypha 
contains a number of effective microtubule tracks which 
run continuously to the tip. Each effective microtubule 
is modelled by a 1-D lattice, with the motors plus cargo 
represented by particles which progress along the lattice 
sites. We identify lattice site 1 with the hyphal tip. Par- 
ticles obey hard-core exclusion and move in one direction 
only, toward the tip, without overtaking. When particles 
leave lattice site 1, they act to extend the lattice through 
the transformation: particle new site. 

We justify our model with a simple order of magni- 
tude test. We identify the lattice repeat with the kinesin 
step size, 8 nm [23| . In the model organism A'^. crassa, 
a 10 fim wide hypha growing at 25° C has as extension 
rate of ^ 0.5 fim per second. It has been estimated that 
~ 600 vesicles per second must fuse with the tip to pro- 



2 

1 2 3 4 .... N 

FIG. 2: Schematic of the model with input rate a, hopping 
rate p = 1, absorption rate /3 and growth rate 7. 



vide enough plasma membrane and other materials to 
maintain this hyphal extension rate 24| . Thus, the ex- 
tension of a hypha by one model lattice unit (~ 10 nm) 
is equivalent to the arrival of order 10 vesicles. If each 
particle delivered to the end contributes to lattice exten- 
sion, we require ~ 10 equivalent effective microtubules in 
a typical hyphal cross section. This is acceptably within 
an order of magnitude of the number of microtubules ob- 
served near the Spitzcnkoi^per within the growing hyphal 
tip of A^. crassa I25j . 

The model dynamics are specified by the rates at which 
the following processes occur on the lattice: particles in 
the bulk hop toward the tip with rate 1; particles enter 
the lattice far from the tip with rate a; particles detach 
from site 1 with rate f3 and transform into a new lattice 
site with rate 7, as shown schematically in Fig. [2l Thus 
7 is the parameter controlling the lattice growth and /3 
allows particles to leave the end of the lattice without 
extending it. Ratio 7//? is thus the efficiency with which 
the hypha extends. Biologically, 7 represents the rate at 
which secretory vesicles fuse with the hyphal tip, a rep- 
resents the vesicle density far from the tip and P allows 
vesicles to reach the tip without contributing directly to 
growth. 

We initially perform Monte Carlo (MC) model simu- 
lations by stochastically updating particles on a single 
lattice according to the above dynamics. After some re- 
laxation time, density profiles are obtained by averaging 
site occupancies over many updates. We find three dif- 
ferent macroscopic behaviours. Results for representative 
parameter values a — 0.25, f3 — and 7 in the range 0.2 
to 0.56, are shown in Fig. [H For high values of 7 one 
sees profiles that decay from the tip to a 7-independent 
bulk density equal to a. For the highest values of 7 the 
density at the tip is < a. As 7 is lowered the tip den- 
sity is > a and the region over which the decay occurs 
grows in size. For low values of 7 we see distinct pro- 
files where the bulk density is 7-dependent and is > 0.4 
(these profiles were also seen for low 7, high a). The 
transition from the high 7 to the low 7 profiles is dis- 
continuous and involves a jump in the bulk density. In 
the regime of high a and 7 (not shown) , density profiles 
with algebraic decays between the boundaries were ob- 
served. These are "maximal current" profiles, which we 
shall discuss shortly. 
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FIG. 3: Average site occupancy for a — 0.25. Upper two 
traces are for 7 — 0.2 (higher) and 0.24 (lower). Lower five 
traces (highest to lowest) for 7 — 0.28 — 0.56 display apical 
peaks and decays. 



We now seek an analytical understanding of our ob- 
servations using a mean-field approximation where we 
consider the average density, pi(t) at site i, and ignore 
correlations between the density at different sites [l9| . 
We describe the growth dynamics in a frame of reference 
co-moving with the growing tip. The tip site is always 
labelled site I. Each time growth occurs, all other site 
labels must therefore be updated i i + The change 
in density at site i is the net result of particles entering 
from the site on the right, departing to the site on the left 
and shifting right due to index relabelling during growth. 
Within the mean- field approximation, we have, for i > 2, 



dpi 

— = pi+i[l - p,,\ - Pi[l- p^ 



-i]+lPi[Pi-i - Pi\ ■ (1) 



Note that the third term is proportional to the rate at 
which lattice sites are added, jpi = v which is the tip 
velocity. Separate equations govern the change in density 
at sites 1 and 2, in order to take into account the effect 
of the particle — > new site transition: 



dpi 
dt 

dp2 
dt 



P2[l-Pi]-(7 + /3)Pi , (2) 

P3[l - P2] - P2[l - Pi] - 7P1P2 ■ (3) 



Eq. [3]differs from the bulk equation only in that should 
a growth event occur, there is never particle at site 2 after 
the lattice indices are updated. The final term is hence 
a decrease in density at site 2. 

Finally, since particles enter at rate a, the particle den- 
sity at the right-hand end is effectively a. As the lattice 
grows, this boundary recedes from the tip with velocity 
—V, ultimately corresponding to boundary condition 



lim pN = a 



(4) 



We seek a steady state solution for this system defined 
in the reference frame of the tip. Such a solution is char- 
acterised by a constant current of particles everywhere 
through the system, a uniform tip velocity and a den- 
sity profile which decays to the right boundary condition 
over a finite length scale, so that the profile is effectively 
independent of the system size, i.e. we seek solutions to 
dlHS]) with the time derivatives set to zero and obeying 
the boundary condition We obtain an expression 

for the particle current through the system in the tip's 
stationary frame from Eq. [1] 



J = Pi[i - Pi-i] - ypi 



(5) 



Now, at the tip J = (7 -|- (3)pi = (1 + l3/^)v, whereas 
far away from the tip ([5]) yields J = q;(1 — a — u),so that 
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(6) 



which gives the tip velocity in terms of a and /3/7. 

We now restrict ourselves to = 0, and comment on 
the effects of non-zero /3 later. Since for /? = the tip 
velocity is simply a result of a flux of particles through 
the final lattice site, we have J — v and thus from ([5]) a 
recurrence relation between the steady state density at 
any site and that at the previous site: 



v{l+P;-l) 
1 - Pi-1 



i > 2 



(7) 



We define poo = a as the stable fixed point value to which 
this recurrence relation converges: 



1 



\/l — 6u + 



(8) 



The decay length to a is finite and independent of lattice 
size, as required. We are now able to solve for all densities 
in terms of the parameters, a and 7: 



Pi 



P2 



7Q! H- 7 

7a [1 — a] 



1 — Pi 7[a + 1] ~ ol[1 — a] 



(9) 
(10) 



and for i > 2, p.; is given through ([7|). 

Under our constraint that the bulk density is a, we 
find two types of steady-state solution to the mean-field 
equations. In these solutions the profiles decay exponen- 
tially toward p = a and are distinguished by whether p 
decays to a from above or below. For 7 > (1 — a)/2, 
which we refer to as region I, the decay is from above 
and for 7 < (1 — a)/2, which we refer to as region II, 
the decay is from below, from a minimum value at site 
2, although there is a peak in the density at site 1. 

These steady state solutions only exist in certain pa- 
rameter regimes. For 7 < a/(l -I- a), instead of iterating 
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FIG. 4: Phase diagram given by the simple mean-field theory 
presented in the text, for the model with /3 = 0. Phases are 
discussed in the text. Sample points for MC simulations of 
Fig. [3] are marked with x's. 

to the fixed point ((H]), the density is fixed for i > f at 
Pi = 1 — 27, with pi = vj^ and v = 7(1 — 27)/(l— 7). The 
interpretation is that the rate of release of particles at the 
growing end is no longer large enough to control the input 
rate. Thus the particle density reaches a maximum value 
that extends from near the tip throughout the whole lat- 
tice and the boundary condition (jl]) is not satisfied. This 
is not a steady state solution for our model in the sense we 
have defined and we describe this region as a "jammed" 
phase. At the transition to the jammed phase the bulk 
density jumps discontinuously from p = a = 7/(1 — 7) 
to p 1 - 27. 

We see from ([8]) that the maximum value of a is q;c = 
\/2 — 1, which is obtained when v = 3 — 2\/2- For a > ac, 
Eq. [T]no longer has real fixed points, and again we do not 
satisfy the boundary condition We may understand 
the region bounded by a > — 1 and 7 > a/{l + a) 
as a maximal current phase, where the particles have 
reached a maximum flow rate through the system which 
is no longer limited by the input and growth rates. In 
this case, the density profile decays algebraically from the 
boundary sites 1 and TV to a bulk density p = ^/2— 1 and 
does not constitute a steady state in our sense since the 
densities evolve as the system grows. 

We summarise the results of this mean-field theory 
with a phase diagram in Fig. [H Regions 1 and II cor- 
respond to the MC profiles observed for high 7 (Fig. [3]), 
and the jammed region corresponds to the profiles ob- 
served for low 7. Simulations carried out over the whole 
parameter space revealed that the qualitative behaviour 
of the mean-field theory is correct, however, the transi- 
tions between different profile types do not occur exactly 
at the predicted mean-field boundaries. In order to com- 
pare in more detail the mean field and MC results we 
plot in Fig. O mean- field and simulation profiles in re- 



FIG. 5: Mean-field density profiles compared to MC density 
profiles for /3 = and /3 > 0. q = 0.25, 7 = 0.3 



gion I. The decay length at the tip is significantly higher 
in the simulation, by a factor of about 10. The differ- 
ences between the mean-field theory presented here, and 
simulations can be attributed to density fluctuations and 
correlations in the system which are ignored in the mean- 
field theory. In particular the occupations of sites 1 and 2 
are strongly correlated since a growth event vacates both 
of these sites simultaneously. An improved mean-field 
theory which takes into account the correlation between 
sites 1 and 2 predicts a phase diagram with essentially 
the same phases, but some modified phase boundaries 




Simulation and mean-field results for non-zero 
Fig. [5l show that /? does not affect the qualitative profile 
shape. A detailed analysis of the phase structure with 
(3^0 will be given elsewhere [28] . 

We take particular interest in the steady state phase 
with a positive density gradient at the tip, and speculate 
whether the emergence of this high density region may 
be associated with the vesicle accumulation within the 
Spitzenkorper in the hyphal tip region of a real fungus. 
For representative parameter values a — 0.25 and 7 = 
0.24, MC simulations predict a high density region of 
^ 240 nm, correlating (within an order of magnitude) 
with the length scale of the Spitzenkoi^per in N. crassa of 
< 2^m [i[27|. 

Note that, in contrast to the vesicle supply centre 
model 0, 0], a growth velocity arises naturally for our 
model. At least within mean-field theory, this velocity is 
determined in a simple way by the rate a at which mate- 
rial is fed into the system far away from the tip and the 
ratio (3/ J representing the efficiency with which the car- 
goes fuse with the tip. It would be of interest to test this 
prediction experimentally, for example by live-cell imag- 
ing the dynamics of motor proteins or cargo fused with 
the green flourescent protein pij . In this way the rate 
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of delivery of motor proteins (or cargo) to the hyphal tip 
could be analysed in relation to the rate of hyphal exten- 
sion (which can be controlled, for example, by changing 
the temperature of incubation) . Finally we mention that 
the phase diagram Fig. S] is related to that of the open 
boundary TASEP |20[|: we explore this correspondence in 
a future publication [28| . 

We thank Graham Wright for many helpful discussions 
on fungal biology. 

KS is funded by the EPSRC. 
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